Bootstrap Percolation on Complex Networks 
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We consider bootstrap percolation on uncorrelated complex networks. We obtain the phase dia- 
gram for this process with respect to two parameters: /, the fraction of vertices initially activated, 
and p, the fraction of undamaged vertices in the graph. We observe two transitions: the giant active 
component appears continuously at a first threshold. There may also be a second, discontinuous, 
hybrid transition at a higher threshold. Avalanches of activations increase in size as this second 
critical point is approached, finally diverging at this threshold. We describe the existence of a special 
critical point at which this second transition first appears. In networks with degree distributions 
whose second moment diverges (but whose first moment does not), we find a qualitatively different 
behavior. In this case the giant active component appears for any / > and p > 0, and the dis- 
continuous transition is absent. This means that the giant active component is robust to damage, 
and also is very easily activated. We also formulate a generalized bootstrap process in which each 
vertex can have an arbitrary threshold. 

PACS numbers: 64.60. aq, 05.10.-a, 64.60.ah, 05.70.Fh 



I. INTRODUCTION 

Bootstrap percolation serves as a useful model to de- 
scribe in detail or in analogy a growing list of complex 
phenomena, including neuronal activity [J-Qjiamming 
and rigidity transitions and glassy dynamics [J, Q , and 
magnetic systems . Chalupa et al. Q introduced boot- 
strap percolation in a particular cellular automaton used 
to study some magnetic systems (for other applications 
see Ref. H), see also the even earlier work of PoUak and 
Riess Q . The standard bootstrap percolation process on 
a lattice operates as follows: sites are either active or 
inactive. Each site is initially active with a given prob- 
ability /. Sites become active if k nearest neighbors are 
active (with k = 2,3, ...). In the final state of the process, 
the fraction 5a(/) of all sites are active. Remarkably, the 
function iSa(/) may be discontinuous. It may have a jump 
at a bootstrap percolation threshold fc2- We will see be- 
low that when this process takes place on a network, this 
is not the only threshold in this system. 

Bootstrap percolation has been thoroughly studied on 
two and three dimensional lattices (see |9l-[l2| and ref- 
erences therein). The existence of a sharp metastability 
threshold for bootstrap percolation in two-dimensional 
lattices was proved by Holroyd M, and later generalized 
to c?-dimensional lattices [10, llli . More recently, boot- 
strap percolation has been studied on the random regu- 
lar graph [l^, [l3| , and also on infinite trees [l^ . Finite 
random graphs have also been studied [ll]. Watts pro- 
posed a model of opinions in social networks in which 
the thresholds at each vertex is a certain fraction of the 
neighbors, rather than an absolute number |l7|. Boot- 
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strap percolation is closely related to another well known 
problem in graph theory, that of the fc-core of random 
graphs |18l421|. The fc-core of a graph is the maximal 
subgraph for which all vertices have at least k neighbors 
within the fc-core. It is important to note the difference 
between the stationary state of bootstrap percolation and 
the fc-core. Bootstrap percolation is an activation process 
which starts from a subset of seed vertices and spreads 
over a network according to the activation rules described 
above. The fc-core of the network can be found as an 
asymptotic structure obtained by a subsequent pruning 
of vertices which have less than fc neighbors. While the 
fc-core has been extensively studied, there are no analyt- 
ical investigations of bootstrap percolation on complex 
networks. 

In this paper we describe bootstrap percolation on an 
arbitrary sparse, undirected, uncorrelated complex net- 
work of infinite size. Specifically, we use the configuration 
model (a random graph with a given degree sequence). 
We show that there are two types of critical phenomena: 
a continuous transition corresponding to the appearance 
of the giant active component, and a second, discontinu- 
ous, hybrid phase transition combining a jump and a sin- 
gularity. (This transition is also often called "mixed".) 
We show that network inhomogeneity strongly influences 
the critical behavior at the appearance of the giant ac- 
tive component in networks with divergent second and 
third moments and finite first moment of the degree dis- 
tribution. In contrast, the hybrid phase transition has 
the same critical singularities for any network with finite 
second moment of the degree distribution. This second 
transition can be understood by considering the "subcrit- 
ical" clusters of the network, consisting of vertices whose 
number of active neighbors is one less than the thresh- 
old. We show that these subcritical clusters give rise 
to avalanches of activations which become increasingly 
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large as the threshold is approached. We also describe 
how the behavior changes when the network is damaged. 
The damaging here is the uniformly random removal of 
vertices, so that a fraction p of vertices in a network are 
retained. We give the phase diagram showing the thresh- 
olds with respect to both the extent of damage to the 
network, and to the size of the initial seed group. In par- 
ticular, there is a special critical threshold Ps at which the 
discontinuous transition first appears. We also show that 
network topology can have a dramatic effect, as on so 
called scale-free networks with finite mean but divergent 
second moment of the degree distribution a qualitatively 
different behavior occurs. There is no phase transition 
in the p-f plane, but instead a giant active component 
appears at any / > and is robust to any amount of dam- 
age (p > 0). Finally, we generalize bootstrap percolation 
by considering a distribution of threshold values, so that 
each vertex may have its own threshold value. We briefly 
outline the equations for the active fraction of the net- 
work and the size of the giant active component in this 
general formulation, and show how the classical percola- 
tion problem and the usual bootstrap percolation (which 
we analyze in the remainder of this paper) are limiting 
cases. 



II. RESULTS 

Consider an arbitrary, sparse, uncorrelated complex 
random network in the infinite size limit. The struc- 
ture of this network is completely determined by its de- 
gree distribution P{q). An important (and convenient 
for analytical treatment) feature of this architecture is 
local tree-likeness, which means that finite loops can be 
neglected. For the time being we assume there are no 
vertices with degree zero. We denote by (q) the mean of 
the degree distribution and similarly (g^) is the second 
moment. The network may also be damaged by the uni- 
formly random removal of vertices so that a fraction p of 
all original vertices remaining. 

Vertices have either an "active" or an "inactive" state. 
Once activated, a vertex remains active. With probabil- 
ity / each vertex is part of the seed group, and is in an 
active state from the start. The remaining vertices, (a 
fraction 1 — /) become active only if they have at least 
k active neighbors. We iteratively activate vertices that 
meet this criterion until a steady state is reached. 



We define to be the fraction of the vertices in the 
graph which are active at equilibrium (that is, includ- 
ing all active vertices, even those forming finite clusters), 
which is also the probability that an arbitrarily selected 
vertex is active in the final state of the bootstrap perco- 
lation process, and the size of the giant active component 
to be iSgc, equal to the probability that an arbitrarily se- 
lected vertex belongs to the giant active component. By 
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FIG. 1: Phase diagram of bootstrap percolation in the f-p 
plane for networks with finite second moment of the degree 
distribution, for k > 2 and smaller than an upper limit femax 
determined by the degree distribution. The solid line marks 
/ci, the continuous appearance of the giant active component 
from 0. The giant component of active vertices is present 
above this line in the upper-right portion of the diagram (la- 
belled II) and absent in the area to the lower-left (I). The 
dashed heavy curve represents the discontinuous transition, 
fc2- This line ends at the special critical point ps- Thin hor- 
izontal lines show the location of the traces in Fig. [2] relative 
to the phase diagram features. 



giant active component we mean a subgraph of active ver- 
tices which forms a connected component that occupies 
a finite fraction of the network. 

In Fig. [T] we show a representative phase diagram for 
the giant active component in the f-p plane for an un- 
correlated infinite complex network whose degree distri- 
bution has finite second and third moments. Results are 
qualitatively the same for any such network. The giant 
active component is absent in the region labelled I, and 
present in the region labelled II. We see that if the net- 
work is sufficiently damaged so that the proportion of 
remaining vertices is less than a critical threshold Pc , the 
giant active component never appears, for any number of 
seed vertices. Above Pc = {q)/[{q^) — (?)], which is equal 
to the well known percolation threshold (see, for example 
[2^), the giant active component appears at some value 
of /j /ci) for a given value of p. This threshold is marked 
by the solid heavy line in the figure. This threshold is 
above zero for all p > pc- In the limit fc — 00 the bound- 
ary between regions I and II tends to the line pf = Pc, 
as the seed vertices may form a giant component in the 
graph. For large values of p, above a special critical point 
Ps, we discover a second transition in the size of the gi- 
ant active component. For a given p > ps there is a 
threshold fc2 > fci at which the size of the giant ac- 
tive component (and also the active fraction) jumps sud- 
denly. These points are marked by the heavy dashed line 
in Fig. [T] Note that Fontes and Schonmann [ij] noticed 
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on undamaged regular random graphs the two transitions 
we observe in more complex networks. 




FIG. 2: Probability that an arbitrarily chosen vertex is (a) ac- 
tive and (b) in a giant connected component of active vertices 
for an Erdos-Renyi graph of mean degree 5, with = 3. The 
four lines (labelled from 1 to 4) are: p = {0.7, 0.893, 0.93, 1}, 
corresponding to the relative positions in the phase dia- 
gram shown as thin horizontal lines in Fig. [T] namely: 1. 
Pc < p < Ps; 2. p = ps; 3. ps < p < 1 and 4. p = 1. The 
arrows in each plot mark the point of emergence (/ci) of the 
giant active component for p = 1. The small dot on each 
trace marks the point of the hybrid transition. 



Fig. [5] shows the active fraction iSa and the size of the 
giant active component iSgc as a function of / for four 
values of p in an Erdos-Renyi graph (which has a Poisson 
degree distributions in the infinite size limit). We choose 
this network as it is representative of random graphs. For 
comparison, the position of each in the phase diagram is 
marked by a faint solid line in Fig. [T] Line 1 is for a 
value of p before the appearance of the jump. Line 2 is 
exactly at the special critical point Ps at which the jump 



appears. Line 3 is at a p > where there is a jump. 
The location of the jump moves to smaller values of / as 
p increases, but never reaches zero, as is demonstrated 
by line 4, which is at p = 1. The giant active component 
appears continuously and linearly from zero, exactly as it 
does in ordinary percolation |22h24|. It is interesting to 
note that there is no discontinuity in iSa at /d (marked 
by a small arrow), i.e., the threshold is invisible (hidden) 
when observing only the overall activation of the network. 

For k = 1 the jump does not appear, we have only 
the continuous transition. For larger k there is a jump, 
and it appears at larger values of / (for a given p) the 
larger k is. The value of Ps also increases, such that there 
is a finite maximum threshold, fcmax (proportional to 
the mean degree of the network for Erdos-Renyi graphs) 
beyond which the jump no longer appears. That is, 
the dashed line in Fig. [1] which marks the location of 
the jump, moves to the right and towards the top of 
the graph as k increases, finally disappearing completely 
above fcmax- This discontinuous transition has a hybrid 
character, with both a discontinuity and a singularity: 
when approaching /c2 from below the size of the giant 
active component approaches the value at the bottom of 
the jump as the square- root of the distance from fc2'- 



S.if) = 5a(/c2) - a(/e2 - /) 



1/2 



(1) 



where a is a constant (see Section |TTI] below for the origin 
of this equation). Lines 3 and 4 in Fig. [2] illustrate this 
situation. The same result holds with respect to p if we 
were to approach this jump along a line of constant /. 
The size of the giant active component, iSgc, has the same 
critical behavior. 

The height of the jump decreases with decreasing p 
(while /c2 increases slightly), disappearing at the special 
point Ps- (The line labelled 2 in Fig. [2] is at p = p^.) 
At this point the behavior is different, as the size of the 
giant component approaches fc2 now as the cube-root of 
the distance from the threshold (see Section llll| : 



5a = 5a(/c2) - a'(/c2 - /) 



(2) 



where a' is a constant. 



To understand the discontinuous "jump" in the size 
of the active component of the graph, we consider the 
subcritical clusters, a concept related to the corona clus- 
ters which were used to describe a similar transition in 
A:-core percolation [25l - [27| . The subcritical vertices of the 
graph are the vertices whose number of active neighbors 
is precisely one less than the threshold of activation for 
that vertex. An example of a small subcritical cluster is 
illustrated in Fig. [31 

Clusters of subcritical vertices are important because 
of the following quality. The activation of even a sin- 
gle vertex neighboring the subcritical cluster necessarily 
leads to at least one of the members of the cluster now 
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FIG. 3: Left: A subcritical cluster in a network with threshold 
k = 3. Filled black vertices are seed vertices, shaded vertices 
are active vertices while empty vertices are inactive. The 
vertices inside the dashed loop form a subcritical cluster, while 
those connected to them are either inactive or have more than 
the threshold number of active neighbors. If a single vertex 
neighboring a subcritical cluster becomes active (for example, 
if it became a seed vertex) — here marked by a cross — it's 
neighbor inside the subcritical cluster must become active, 
and then it's neighbors, until the entire cluster is activated, 
as shown on the right. 



meeting its activation threshold. In turn, this will ac- 
tivate one of its neighbors in the cluster, and so on, so 
that an avalanche of activations ensues until the entire 
subcritical cluster becomes active — see Fig. |3l Below and 
above the jump, the subcritical vertices form only finite 
and isolated clusters, but as fc2 is approached from be- 
low, the mean size of the subcritical clusters diverges. 
Hence the avalanches resulting from the change in ac- 
tivation state of a single vertex form a finite fraction of 
the entire graph, leading to a discontinuous change in the 
size of the active fraction of the graph. This argument 
will be made more precisely in the analysis in Section IIVI 
below. 

When the degree distribution of the network decays 
very slowly, specifically in networks with divergent sec- 
ond moment (g^) or third moment (g"^) of the degree dis- 
tribution, the results are different from those described 
above. In particular, this is the case if the degree dis- 
tribution tends to the form P{q) oc q~'^ with 7 < 4 for 
large q. If 7 < 2 the mean of the degree distribution 
also diverges, but we do not consider this case here. If 
2 < 7 < 3, the second moment diverges. These scale-free 
networks are of particular interest because many large 
real natural and technological networks appear to be of 
this kind [2^, [IJ]. In this case (in the limit N — s- 00) a 
finite fraction of vertices is activated at any p or / > 0, 
and for any arbitrary activation threshold k. In other 
words the location of the jump tends to zero as the size 
of the network increases — so in very large scale-free net- 
works we will not find (at finite / or p) either of the 
transitions observed in graphs with fast decaying degree 
distributions. This means that, as has been found in sev- 
eral other cases [2^, [2^ , such scale- free networks are very 



robust to damage, and also that such a network is very 
easily activated. 

When 3 < 7 < 4, the phase diagram is qualitatively the 
same as that shown in Fig.[T] The giant active component 
appears at finite p ( or /) with a continuous transition. 
However, rather than growing linearly near the transition 
point, the size of the giant active component increases as 
the distance from the critical point raised to the power 
1/(7-3). 



III. BASIC ANALYSIS 

In this section and the two following, we describe in 
more detail how the results already described may be 
obtained. 

Consider choosing an arbitrary vertex from the net- 
work. We wish to calculate the probability that this 
vertex is active in the equilibrium state. To calculate this 
probability, we first define Z as follows: Z is the proba- 
bility that, on following an arbitrary edge in the graph, 
we reach a vertex which is either a seed vertex or has at 
least k downstream neighbors that are active. (By down- 
stream we mean neighbors of the vertex reached by the 
edges other than the one we arrived from.) To be ac- 
tive, these downstream neighbors in turn must fulfil this 
same condition, that they are either seed vertices or they 
have k further downstream neighbors of their own that 
are previously active. 



• o 

seed vertex non-seed vertex 




O © -^gc 

TABLE I: Symbols used in graphical representations of self- 
consistency equation. 

We can graphically represent this recursive relationship 
using the symbols given in Table HI The probability Z 
is represented by an edge ending in a square. A seed 
vertex is represented by a black disc, and other vertices 
by open discs. An edge crossed by a short line at its end 
represents the probability 1 — Z, that is the probability of 
encountering a vertex that doesn't satisfy the condition 
for Z. Thus we obtain the following representation for 
Z: 
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The terms on the right hand side represent sums of the 
probabiUties of all such terms. Based on this diagram, 
we write mathematical expressions for the probabilities 
represented by each of these symbols, allowing us to con- 
struct the following self-consistency equation for Z: 

^ 

>^jz{\^z\i~zr-K (3) 

The probability iSq, represented by a shaded disc, is the 
sum of two terms, as represented in this diagram: 



The first is the probability that the vertex is active from 
the beginning (p/), the second [with pref actor p(l — /)] is 
the probability that it has at least k neighbors that would 
be active even if the vertex we are observing was inactive. 
But each neighbor satisfies this condition precisely with 
probability Z, as represented by squares in the diagram. 
Converting to a mathematical expression, this gives the 
following equation: 

OO i / .\ 

S,=p}+p{l-})Y^P(^)Y^r\z\l- Zy-' . (4) 

i=k l=k ^ ^ 

The probability 5gc that an arbitrarily chosen vertex 
belongs to the giant active component can be constructed 
in a similar way, but we must impose the further condi- 
tion that a vertex has an edge leading to an active sub- 
tree of infinite extent. We define X to be the probability 
that the vertex encountered upon following an arbitrar- 
ily chosen edge meets the conditions for Z and also has 
an edge leading to an active subtree of infinite extent. 
Graphically, we represent the probability X by an infin- 
ity symbol at the end of an edge, and a self-consistency 
condition for X is expressed by the diagram: 




This corresponds to the equation: 

i=0 m=l ^ ^ 

X ^ (5) 

m=l ^"^^ 

The probability Sgc (the probability that an arbitrary 
vertex belongs to a giant active component, represented 
below by a shaded circle containing the infinity symbol — 
see Table H]) is the sum of the probability that the vertex 
is a seed vertex that is connected to an infinite active 
subtree (probability X) and of the probability that the 
vertex is not a seed vertex but has at least k indepen- 
dently active neighbors (probability Z), at least one of 
which also leads to an infinite active subtree. Thus 




so that, similarly to Eq. we can write iSgc in terms 
of X and Z: 

OO i 

z^O m— 1 

OO i / 

i=k l=k ^ ^ 

^ E ( ^ ) x"(z - x)'-"(i - zy~K (6) 

m=l V"^/ 

It will be useful to define 'i'{Z,p, f) to be the right-hand 
side of Eq. and ^{X, Z,p, f) the right-hand side of 
Eq. ([5]), so that these two equations become 

Z-v|/(z,p,/), (7) 

and 

X = ^X,Z,pJ). (8) 

These equations can be solved numerically for a given 
network degree distribution. If multiple solutions exist, 
the physical solution for Z is always the smallest value. 
The location, /d of the appearance of the giant active 
component can be found by assuming X is small but 
non-zero in Eq. ([8]), taking the limit as X tends to zero 
and solving for / for a given p (or vice versa). In this 
way we also find that X and hence iSgc grow linearly from 
the critical point /d for networks with (g^) finite. The 
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results mentioned below also correspond to this case — we 
will examine the case (g^) — >■ oo subsequently. 

The second, discontinuous, transition can be located 
by observing that the jump occurs when, (after a second 
solution appears) the smallest Z solution of Eq. ([T]) disap- 
pears. At this point '^{Z) just coincides with the value 
of and a little consideration reveals that this must 
be at a local extremum of jZ. Thus the discontinuous 
transition can be found by simultaneously solving Eq. ([7]) 
and 



dZ\Z 



= 



(9) 



for /. The fact that the first derivative vanishes leads 
to the square-root scaling near the critical point, with 
respect to either / or -p — see Eq. ([!]). The jump dis- 
appears at a special critical point ps in the f ~ p plane 
which satisfies Eqs. (O and a third condition 

dZ^ 



(10) 



This condition means that the scaling below (see 
Fig. [T]) is cube- root — see Eq. 



IV. AVALANCHES 

The singular behavior [Eq. ([ij] near the hybrid tran- 
sition can be understood by considering the subcritical 
clusters of the active subgraph. These form a subset of 
the inactive portion of the graph consisting only of those 
vertices whose number of active neighbors is exactly one 
less than the activation threshold for that vertex — see 
Fig. [3l That is, the subcritical subgraph consists of all 
those vertices which are not seed vertices and which have 
exactly A; — 1 active neighbors. The subcritical clusters 
are finite everywhere except exactly at the point of the 
jump transition. To show that this is the case, we use 
agenerating function approach, similar to that used in 
HE US to calculate (s)sub, the mean size of the sub- 
critical clusters. 

Let ^0(2;) be the generating function for the proba- 
bility that an arbitrarily chosen vertex is subcritical. By 
considering the probability that an arbitrarily chosen ver- 
tex is subcritical, which corresponds to Fo{l), we can 
write 



Foix) =p{l-f) 



q>k-l 



J2 Piq)L^^^Z''-\l-Z)'^-''+'x'^-''+K 



(11) 

Similarly, the generating function for the probability that 
an arbitrarily chosen edge leads to a subcritical vertex is 

F, (.) = Ki - /) E TT (r 1) zy-'-'-'- 



q>k 



(12) 



The generating function for the probability that a ran- 
domly chosen vertex belongs to a subcritical cluster of a 
given size then must obey the self-consistency equation 

Mini 



Ho{x) ^ 1 ^ Foil) + xFo[Hi(a 



(13) 



where 1 — Fo{l) represents the probability that the ran- 
domly chosen vertex is not itself subcritical, and the sec- 
ond term is a recursive relationship, ensuring that suc- 
cessive powers of x correspond to the probabilities of en- 
countering a cluster size matching that power. In this 
equation Hi (x) is the related generating function for the 
probability that a subcritical cluster of a given size is 
reached upon following an arbitrarily chosen edge. In a 
similar way we can write a self-consistency equation for 
this: 



Hiix) ^ 1 - Fi{l) + xFi[Hiix)]. 



(14) 



Where 1 — i^i(l) is the probability that the arbitrarily 
chosen edge leads to a vertex that is not subcritical. Note 
that Ho{l) = Hi{l) = 1. 

From these generating functions, we can calculate var- 
ious quantities related to the subcritical clusters. For 
example, the distribution of avalanche sizes (which are 
the same as the sizes of the subcritical clusters) is given 

by 



1 d^Hojx) 
G{s) = 



(15) 



x=0 



and we expect that at the critical point G{s) ~ s 
The mean size of the subcritical clusters is simply 



dHn 



dx 



x = l 



q>k-l 
7k — lf-\ v"!?— 



-^Fo{l)+p{l-f) E P^9\k_l 

dHi 



X Z^-'il - Zf 



(g-fc + l)- 



dx 



(16) 



x = \ 



Using Eqs. and comparing with Eqs. and 

© we find that 



dHx 



dx 



\-d'^{Z)/dZ' 



(17) 



Now from Eq. ([9]), d^{Z)/dZ = 1 at the critical point, 
and 1 - d-^{Z)/dZ (X Z - Zj oc (/c2 - fY^^, near the 
critical point. Thus, near this point, therefore, the term 
containing (i_ffi/c?x|^ .|^ dominates (s)sub so that 



(s)sub « (/c2 - /)-l/^ 



(18) 



or alternatively, for fixed /, (s)sub {pc2 ~ p)^^^^, hence 
the mean size of the corona clusters diverges at the crit- 
ical point. 
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The addition of a single vertex (an infinitesimal in- 
crease in p) or activation of a seed vertex (increment 
of /) may lead to the activation of a subcritical ver- 
tex and hence activating an entire subcritical cluster in 
an avalanche. At fc2, the subcritical clusters span the 
whole graph, so the activation of a vertex can lead to 
an avalanche of activation that eventually affects a finite 
fraction of the whole infinite graph hence we see a dis- 
continuity in both the size of the active fraction and the 
giant active component. Note that for / > /c2 the mean 
size of the subcritical clusters is finite. 



V. SCALE-FREE GRAPHS 

Let us consider degree distributions that tend to 
P{q) cx for large q where exponent 7 is some positive 
constant usually > 2. For concreteness, in the following 
we will consider the degree distribution 



Piq) = Aq- 



for q>qo, 



(19) 



where ^ is a constant of normalization. For 7 > 4 the re- 
sults are qualitatively the same as those described above. 

When 7 < 3 the second moment of the degree distribu- 
tion diverges, leading to different behavior. The results 
that follow refer to the situation when 2 < 7 < 3. Note 
that many real world networks, especially biological net- 
works have exponent in the range 1 < 7 < 2 [23^ [23, H^] . 
In this case the first moment also diverges. We don't 
consider this case here. 

By assuming Z to be small, we can approximate '^{Z) 
by considering only leading order in Z. Then the self 
consistency equation ^ becomes: 



Z^pil-f)aZ^-^+pf. 



(20) 



For p > 0, this equation has no small-Z solution, even 
in the limit / — >■ 0. Because the only solutions for Z 
as / — 7- are therefore of order 1 it is clear that the 
discontinuous transition is moved to / = for scale- 
free graphs. A similar analysis for the giant active 
component — approximating $(Z, A) [the RHS of Eq. 
([5])] by assuming X and Z both small leads to similar 
conclusion about X: that there are no infinitesimal so- 
lutions in the limit / — > 0, confirming that there is no 
jump for / > but also that the giant active component 
appears for any / > and p > 0. 

To add support to this approximation, consider the 
same degree distribution as before, but truncated at some 
maximum degree (/cut (the normalization constant will 
also necessarily change). If we re-derive (|20p assuming a 
finite (7cut, we find 



Z^p{l-f)bZ'^+pf 



(21) 



which does have a solution at finite p (or /). For finite 
gcut, numerical solution of Eq. shows a jump appears 
at small values of /. As (/cut is increased, the curve of 



this jump moves closer to / = 0, and extends towards 
p = 0. Similarly the giant component appears at smaller 
and smaller values of p and / as (/cut is increased. In 
keeping with the approximate analysis just described, we 
expect that both thresholds reach f — and p = in 
the infinite size limit. In summary, when 2 < 7 < 3, the 
giant active component is always present everywhere in 
the p — f plane for p > and / > 0, and appears not 
from zero but with a finite size. 

When 3 < 7 < 4, an expansion of the right-hand-side 
of Eq. ([5]) in leading powers of X gives an equation of the 
form 



X ^ciX + C2A ' 



(22) 



where the ellipsis signifies further terms of higher order in 
X. The first coefficient ci = pf^^^}^ +p{l - f)B{Z), 
where B{Z) is a function of the variable Z . Thus when 
/ < 1 the value of ci differs from that found in the per- 
colation problem. The presence of a finite linear term 
(ci > 0) means that the appearance of the giant compo- 
nent occurs at non-zero values of p (or /) - at a point 
which can be found by solving ci = 1. However, because 
the second leading exponent is 7 — 2 and not 2, X scales 
as {p — Pci)^ , near the appearance of the giant compo- 
nent, with /3 = 1/(7 — 3). This is the same scaling as 
found in the usual percolation problem Curiously, 
the second coefficient C2 is simply equal to p up to a fac- 
tor depending on the degree distribution. Above 7 = 4, 
we find (3 — 1 as found for the usual percolation. 



VI. GENERAL DISTRIBUTION OF 
ACTIVATION THRESHOLDS k 

The bootstrap percolation process described above can 
be thought of as a specific case of a more general pro- 
cess in which the threshold values can be different for 
each vertex. Assuming no correlations between vertex 
degree and threshold value, we can define a distribution 
Q{k) such that Q{k) is the fraction of vertices which have 
threshold value k. The fraction of seed vertices is then 
g(0). Setting Q(0) = / and Q{k) = 1 - / for some 
k > 2 we recover the bootstrap percolation model de- 
scribed above. 

In the general case, we find that the equation for the 
active fraction is 



5a=p5]Q(fc)^P(*) 



k>0 



i=k 



.l=k 



(23) 



where, as above, Z is the probability of encountering a 
vertex with at least k downstream active neighbors upon 
following an arbitrary edge: 



z->-EqwE "^^",rr" E(;)z'('-^) 

fe>0 i=k 



l=k 



(24) 
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These two equations are similar to those presented in [3J| 
for undamaged networks as a generahzation of the Watts 
model [13 ■ 

Similarly, the equation for the giant active component 

is 



k>0 



l=k 



(25) 

where as before X is the probability that an edge leads 
to an infinite active subtree: 



X=pY,Q{k)Y. 



(^ + l)P(^ + l) 



k>Q i=k 

i=k ^ ^ ^ ^ 



-I 



(26) 

Vertices which have k = 1 become active if they have 
a single active neighbor. Thus a single seed vertex will 
activate an entire connected cluster of such vertices. In 
particular, if there is a giant connected cluster in the 
network (i.e. Hp > Pc, the percolation threshold), the 
introduction of a finite number of seed vertices into the 
infinite network will (almost surely) activate the giant 
connected component. In other words, we have behavior 
equivalent to ordinary percolation. In particular, if we 
set (3(0) + g(l) = 1 (and requiring that, if Q{0) -J> 0, the 
number of seed vertices remains sufficient to activate the 
giant component of the network) then we recover from 
Eqs . (pSj) and ([^ the well known percolation equations 



VII. DISCUSSION 

In this paper we have extended the understanding of 
bootstrap percolation to uncorrelated infinite random 
graphs with arbitrary degree distribution, and studied 
the effects of damage to the network. We have found that 
the phase diagram for the giant active component with 
respect to damage to the network {l—p) and the fraction 
of initially active vertices (/) has several interesting fea- 
tures. There are two transitions observed. At the first the 
giant active component appears continuously from zero, 
and at the second (always at a higher initial activation 
fraction) there is a hybrid phase transition, where the 
size of the giant active component has a discontinuity — a 
'jump' — while also having a singularity, as the size of the 
giant active component approaches the transition from 



below as the square root of the distance from the criti- 
cal point. This singular behavior is due to avalanches in 
the activation process. The sizes of avalanches of activa- 
tion are determined by the size of subcritical clusters — 
clusters of vertices whose number of active neighbors is 
exactly one less than the activation threshold. Every- 
where but at the hybrid transition these subcritical clus- 
ters are finite (though together occupying a finite fraction 
of the network), but as the transition is approached these 
clusters grow as the reciprocal of the square-root of the 
distance from the transition. We also observe a new spe- 
cial critical point, at the level of damage at which the sec- 
ond transition first appears. Here the height of the jump 
tends to zero, and the scaling near the critical point is the 
cube-root of the distance from the threshold. These re- 
sults are valid for arbitrary degree distributions, so long 
as they decay rapidly enough that the second and third 
moments of the distribution are bounded. Note that we 
could express our results not in terms of / and p, but of 
/ and any other convenient parameter, for example, the 
mean degree (q) of a network. This allows one to apply 
our conclusions to arbitrary uncorrelated networks. 

Network inhomogeneity plays an important role. 
When the second moment of the degree distribution is 
bounded but the third moment is unbounded, the criti- 
cal scaling near the appearance of the giant active com- 
ponent is not simply linear but has higher order scaling, 
depending on the degree distribution. When the second 
moment is unbounded, for example in scale-free networks 
with degree distribution exponent 7 < 3, both thresh- 
olds tend to / = and p = in the infinite size limit. 
Thus the phase diagram is featureless, with a giant ac- 
tive component (albeit sometimes very small) present for 
any finite activation and any amount of damage to the 
network. This result has important implications for real 
world networks. For example, the network of neurons 
in the brain may have such a scale-free organization 35 1 
meaning that brain activity may be able to be instigated 
with very small stimulus (even though such networks, 
whilst large, are of course finite). 

In summary, we have obtained phase diagrams for the 
bootstrap percolation problem in a wide range of com- 
plex networks. We have described the properties and the 
nature of two distinct transitions in this problem: the 
bootstrap percolation transition and the emergence of a 
giant connected component (percolative cluster) of active 
vertices. 
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